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FOREWORD 


For the past five years the Aircraft Propulsion Laboratory, Department of 
Aeronautical Engineering, University of Michigan, has made extensive experiments 
in interferometric techniques using a Mach-Zehnder Interferometer to study the flow 
of air through a single stage, radial turbine. These experiments were sponsored by 
the Large Steam Turbine Division of the General Electric Company, Schenectady, 
New York. 

The above experiments resulted in several series of interferometric photo- 
graphs for various running conditions of the turbine. Particularly interesting was 
the repeated shedding of small] vortices from the trailing edges of the nozzle blades, 
at generally regular intervals. It was decided to analyze a series of these pictures, 
in an attempt to determine the nature of these vortices, and the accompanying flow 


patterns. This analysis, and the resulting conclusions, are the work of this report. 





INTRODUCTION 


The analysis of the flow in a single stage radial turbine by interferometric 
methods was undertaken in an attempt to determine the frequency and strength of 
small vortices, shed off the trailing edge of the nozzle blade. The flow was considered 
isentropic and the analysis made use of a sequence of 24 flash pictures taken for con- 
secutive turbine blade positions. Thus it was assumed that the flow repeats itself for 
each revolution of the turbine wheel, if the total pressure and flow rate are constant. 
The pictures used, Figs. 1 to 24, were part of a series included in Ref. 2, which was 
generally a feasibility study for the interferometric technique with regard to a moving 
turbine. 

The analysis was carried out during the academic year 1959 - 1960 at the 
Propulsion Laboratory, NOE Campus, University of Michigan, Ann Arbor. The 
work was sponsored jointly by the General Electric Company of Schenectady, New York 
and the United States Naval Postgraduate School, Monterey, California. Lt. L. G. 
Hangartner and Lt. D. H. McVay were the analysts although Mr. H. Kraft of General 
Electric Company and R. E. Cullen, of the University of Michigan, were close ad- 


visors in the project. 





ANALYSIS 


I GENERAL 

In obtaining the interferometric photographs used for the analysis in this re- 
port, an 8-inch Mach-Zehnder Interferometer was employed. The turbine was a radial 
flow, single-stage impulse turbine designed and fabricated by the General Electric 
Company. This turbine has four cut-out positions, every 90°, located in the end 
plates. Thus, the air flow through these four bucket channels, is visible and may be 
photographed. Ref. 1, presents the description of this interferometer, turbine and 
accompanying apparatus. Figures from Ref. 1, and a brief discussion of the principles 
of interferometry may be found in Appendix I to this report. 

Appendix I, also points out the difference between fringe displacement and in- 
finite fringe type pictures. Although both were used in this report, the infinite fringe 
type picture was used much more extensively. Fig. l and 27, are examples of the 
infinite and fringe - displacement pictures, respectively. 

These pictures were taken using an achromatic focusing lens and a 4 in. X 5 in. 
camera back, together with an air cooled mercury arc eraieen apparatus. It was 
also possible to take pictures of the turbine in continuous operation, using a Fastax 
type of 16 mm. motion picture camera. 

Figs. 1 to 24, corresponding to UM 435 to 458, respectively, are the inter- 
ferometric photographs used in the analysis of this report. Assuming that the flow in 
the turbine repeats itself every revolution, it was possible to use these 24 flash pic- 


tures, taken over a period of minutes, but under the same turbine operating conditions, 





and establish a true time study of a turbine blade passage. The flash apparatus was 
electronically synchronized, as discussed in Ref. 1, to "stop' the turbine, as the cut- 
out position progressed across the test section area. Using turbine RPM, held con- 
stant, in this case at 900 RPM, it was then possible to eonvert a reference distance to 
true time. 

As mentioned previously the flow was considered isentropic. The air flowing 
through the turbine is assumed to be a perfect fluid. This assumption of a perfect 
fluid eliminates the property of viscosity found in actual fluids, due to internal friction, 
and the perfect fluid, therefore, can exert no shearing stress. The assumption of a 
perfect fluid also implies irrotational motion in the fluid. In addition, the flow was 
considered to be incompressible throughout the nozzle passage. 

Due to the assumptions above, the standard equation of state and Bernoulli 
equation may be used to compute the velocities around the blade from the density 
measurements of the interferogram. As a result of the irrotational property of the 
flow, a velocity potential may be found, and using this potential in conjunction with 
Kelvins' equation, makes it possible to determine the final velocity profile about the 


nozzle blade. 








II DETERMINATION OF VELOCITIES ON THE NOZZLE BLADE 


In analyzing the pictures, Figs. 1 to 24, the following data taken from Ref. 2, 


was used: 
Gap setting .133 inches 
RPM 900 RPM 
oe 2on02 Ain, Hy abs.) 
ae Ff, Oeste Op 
fo 1.41 (in. HO) 


Converting the above to pressures in psia, and from the equation of state, determin- 


ing the free stream density, the following values were obtained. 
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Using these free~stream state functions, it was then possible to determine the 
free-stream velocity, using the assumption of incompressible flow: 
ee = a So Va / 9 
Ye 2 
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In order to find velocities along the nozzle blade, corresponding to each interfero- 
metric fringe it was necessary to establish the change in density for each fringe shift. 


From Ref. 2, and as explained in Appendix I, this change in density for one fringe- 


shift was given as: 


AS /. _ 
J fe ore. 





For a half fringe shift from free stream conditions, then: 
/ 
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Using Bernoulli's equation, it was possible to solve for the veloeity corresponding to 


each half fringe-shift. 
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In order to find Py> it was necessary to used the isentropie relation: 
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Then the velocity was computed from: 
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Table I, gives the velocities found in this manner for cach 1/2 fringe-shift up to seven 
complete fringe-shifts. 

It was now possible to determine the aetual velocities around the nozzle blade. 
Referring again to Figs. 1 to 24, the blade on the right hand side was chosen as the 
one to be analyzed. Measurements around the blade were made from the trailing edge 
in a countereloekwise direction. The white area in the upper part of the Figure is the 
free-stream air flow. Beginning at the leading edge, the first black fringe indicates 
a half fringe-shift. Each consecutive black fringe indicates a full fringe-shift. Vel- 
ocities corresponding to these fringe-shifts, as found in Table 1, are then assigned to 
positions on the blade at the mid-point of the black fringe. The figures were enlarged 
about 1.75 times and actual positions on the blade were found using the known distanee 


around the blade of 4.98 inches. Table I] is a sample of this procedure. Table II 





contains the final velocities and their true positions on the blades. 

Table III also includes reference time, © , for each picture. This reference 
time was found using the reference distance from a permanent vertical wire down the 
center of the test section, to the mid-chord point on the suction face of the turbine 
bucket. This distance was converted to actual distance, and then converted to true 
time, uSing turbine RPM and the radius to the center of the turbine bucket. 

RPM - 900 
Rad. of turbine wheel-20. 875" 
Chord of turbine blade -. 908" 
Rad. to center of blade-20, 421" 
AD 
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Table IV lists the values of the reference time,* ;for each of the photographs, (where 
Fig. 1 was arbitrarily chosen as U = a), and the computations involved. 

It should be noted in Table III that the leading edge stagnation point is listed at 
the same blade position for each photograph. Before assuming this point, three studies 
were made. Plots of V vs. s were made, of the results in Table III, omitting the 
Stagnation point, and allowing the curve to pass through zero velocity. Examples of 
these curves, together with theoretical curves established by the equations found in 
Ref. 3, are shown in Figs. 25, and 2G. The values of s, where V passed through Zero, 


were found to extend from s = 2.635" tos #« 2.78". The average of the 24 pictures 





however was found to be s = 2.715", 

Secondly, the fringe displacement type interferograms, using the same turbine 
conditions and gap setting, were studied. Fig. 27, is an example of these interfero- 
grams and the series used included UM 413 to 434. Fig. 27 is photo number UM 413. 
A fringe in this case, which bends away from the vertical, indicates a finite velocity. 
Determining either the vertical fringe, or the point between two very slightly displaced 
fringes, indicated the location of the stagnation point. Table IX lists the position of 
the stagnation point as found in both of the above methods. The average of the fringe 
displacement pictures was s = 2.741"'. The average of the two methods was 2.728" 
taken as 2.73''. It was this value which was used in Table II. 

The third method was somewhat a check of the other two. In reviewing Figs. 

1 ~ 24, it was noted that in several cases, a faint fringe may be seen on the leading 
edge of the nozzle blade. By taking the geometrical center of this fringe, a fairly 
accurate position of the stagnation point may be located. Figs. 8, 9, 15 and 16 are 
examples of this, and in the pictures checked, a good corroboration to the value 


s = 2.73" was obtained. 








III Non-Steady Analysis 
When the flow field around the nozzle blade was originally analyzed, the Ber- 


noulli Equation, 
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was used to determine the velocity at any point on the blade. However, in non-steady 
flow of the type encountered in the turbine, it is necessary to use Kelvin's Equation as 
found in Ref. 4: 
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where cis the non steady potential term. To explain this term a little of the theory 
of fluid flow will be presented. The term field denotes a region throughout which a 
quantity is defined as a function of location within the region and time. If the quantity 
is independent of time, the field is steady or stationary. 

The field of gravitational force is a conservative field, meaning by this that the 
work done by the weight in taking a body from point A to another point B is independent 
of the path taken and depends on the vertical height between A and B only. A conser- 
vative field of force gives rise to potential cnergy, which is measured by the work 
done in taking a body from one standard position to another, independent of the path. 


The velocity in a field is assumed to be the result of just such a potential field. 


While this is not physically true, it is an exact mathematical analogy. We may write: 
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where the function @ is called the velocity potential. The velocity potential is, in 


general, a function of (x,y,t), and may be determined from: 





provided the velocity is known. Since the velocity potential is a function of time, it 
necessarily follows that the velocity is a function of time. 

To determine the strength of the shed vortices, the exact velocities in the 
field must be known. To determine the actual velocity, it was then necessary to de- 
termine the magnitude of the CLS term. 

For this analysis, only part of the series of pictures was considered. Figs. 4 
through 20, were chosen, Since they covered two turbine blade passages, and showed 
a clear case of the vortex shedding phenomenon. Since there was some question as to 
the probability of determining the exact magnitude of the ae term, a feasibility 
study was made. For several reasons, in this study, only the suction side of the 
nozzle blade was used. One reason was lack of information on the pressure side due 
to the very few fringes distinguishable. Secondly, the time and effort involved, even 
for a very limited analysis, was considerable when all the calculation is done manually. 
And thirdly, the experimental curve closely approximated the theoretical curve along 
the suction side, leading the investigators to believe, that limiting the investigation 
would improve the results. 

Figs. 28 to 43, are plots of velocity versus perimeter distance, for the suction 
side of the nozzle from s = .25 in. tos = 2.73 in., the assumed stagnation point, for 
the series of pictures to be analyzed. The area from s = 0 tos = .25'' was not con- 
sidered due to the problems encountered not only because of the wake and shed vortices 
in this area, but also because of the build-up of the boundary layer in this region. 

Assuming a potential exists (irrotational flow), the line integral of the velocity 


component in the direction of the path between any two points in an irrotational field 





is equal to the difference in potential between these points. That is: 


re) 
fo Vd = o.- 4, 


A 
If point A is a position of zero potential, the line integral from A to any point in the 


field is the value of the potential at that point. 

Beginning where V = 0, at the nose stagnation point, or at s - 2.73 in., these 
figures were graphically integrated to obtain the potential function @. This integration 
was made every S = .25 in., and was additive along the periphery of the blade. Values 
of and V versus s are tabulated in Table V. These values of 9 and V were plotted 
versus time, using the reference time for each of the pictures, for each of the even 
intervals along the periphery, i.e. .25 inches, .50 inches, etc. These plots are 
found in Figs. 44 to 53. 

Taking the slope of the @ versus t curve at any particular time, gave the in- 
stantaneous value of °?/y¢. Using this value of °?/2t in the Kelvin Equation with an 
iteration process, it should be possible to determine the actual velocity distribution 
about the blade. 


The iteration is based on the general formula: 
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Where the F(t) of the Kelvin Equation is assumed to be a constant, as suggested by 
Lamb in Ref. 4. 
The zeroth iteration or initial determination of the velocity is based on the 


premise Aas ~ 0, which reduces the equation to the Bernoulli equation. Hence, 
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Where Vo is the velocity determined from the interferometric pictures and tabulated 


in Table IU. 


For the first iteration, 


oO”, + mle + Me = ae 
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For this iteration Cj is assumed to equal to C, and throughout the analysis p/y 


is assumed to be constant. Therefore, 
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Since trial calculations showed that the term 2 Sas. was quite large and tended 
to make the values of V oscillate widely, hence making the iterative process extremely 


A 
lengthy, it was decided to use one-half this value or just ~ ise to speed convergence. 
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For the second iteration, 
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It may be shown mathematically that the constants Cg. C,, Co, etc. are all 


equal, and henceforth, will be known simply as C. Therefore, 
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Cancelling Pie and solving for Vo, 


likewise it can be shown: 








or in gencral terms: 
| eo 


Vn *#/Vo - So 

This iteration process has been carried out several times,in an attempt to 
determine the rate of convergence and the approximate magnitude of the final velocity. 

The original velocities, as stated before, are tabulated in Table II and the a) 
and V versus time plots are Figs. 44 - 53. For the first iteration, the original vel- 
ocities, the value Bo determined from the @ plots, and the corrected velocities 
are tabulated in Table VI, 

From these corrected velocities new plots of velocity versus perimeter distance 
were made and are included in the report as Figs. 54 to 61. These figures were again 
graphically integrated and the the values of Pand V versus s tabulated in Table VII. 
From this data, plots were made of the new Dp and V versus time and are included in 
the report as Figs. 62 to 71. From these plots the new values of oe were found 
and the corrected velocity again computed. The tabulated values of d t/ ot and the 
corrected velocities are included as Table VII. 

At this point three iterations had been completed and the three time histories 
of velocity for a specified peripherial distance (s) are plotted in Figs. 72 to 74. 

During the iterative process the question of sign for the MY y term arose, and 


for all the calculations the term was taken as plus. That is, negative slopes on the D 


ata E 





versus t plot were - Obyy, . 

This has by no means been a complete analysis of the non~steady flow in the 
turbine, for many reasons, foremost of which is the limited time available. But the 
work is considered very useful and the foundation for more advanced and continued 
work has been laid. The problems encountered, decisions made, conclusions drawn, 
and recommendations made will all be covered in detail in the discussion and conclu- 


sions sections of this report. 
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DISCUSSION 


The preceding analysis outlines the general method used to reduce the inter- 
ferometric data obtained from a picture series, for a given inlet nozzle blade in the 
turbine. It was hoped that the analysis would provide a clear quantitative indication of 
the strength and frequency of the vortex, which is periodically shed from the trailing 
edge of the blade. This vortex may be seen in Fig. 1, just below and to the left of the 
trailing edge, of the right-hand blade. 

It was the general intent of the complete analysis to obtain a plot of circulation 
around the blade versus time. By analyzing the changes in circulation, which were 
believed to be affected by the vortex shedding, it was thought that a good indication of 
the strength of the shed vortex could be found. Also it was believed that the circulation 
versus time plot would be sinusoidal in nature, and give a good indication of the fre- 
quency and actual times of shedding. 

' From the above data a gencral idea of the power loss associated with the vortex 
shedding phenomenon could be gained, and compared with any experimental data avail- 
able on this subject. 

The results, however, from this limited analysis have not explained, nor sub- 
Stantiated, the pictures of the vortex shedding. Some of the reasons for this deficiency, 
and recommendations for future investigations are the main subjects of this discussion. 

Some idea of the frequency of the vortex shedding was gained from studying 
movies of the flow, using a high speed motion picture camera. It appears that a vor- 


tex is shed from the trailing edge of the nozzle, every time the turbine bucket leading 
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edge passes the nozzle tip. This helped to fix the frequency of the vortex shedding, 
although it was of little help in determining the strength. 

In analyzing this particular sequenee, the blade on the far right of the test sec- 
tion was chosen. The main reason for this choice, was the fact that, this blade was 
mounted most nearly parallel to the flow. It should be noted, the nozzle blades are 
mounted in a slight arc so as to conform to the radial shape of the turbine, thus main- 
taining a eonstant gap setting, between stationary and moving blade rows. The test 
section is just to one side of the center-line of this are. The flow entering the nozzles 
is laminar, and parallel to the referenee wire, mounted in the center of the test sec- 
tion. Henee, the right hand blade was chosen, in order to eliminate any correetions 
necessary, when the blade is mounted at an angle to the flow. 

It is coneeded that there may be advantages to analyzing the flow around other 
blades, partieularly the eenter-most blade, i.e. the second from the right. For the 
series considered, this blade elearly showed the shedding of the vortex, plus giving a 
more detailed picture of the complete flow field, due to its location in the eenter of the 
test section. The error, resulting from the blade being slightly misaligned from the 
direetion of flow, is assumed to be negligible compared with measurement errors. 

In reality, it would be desirable to have a complete analysis of all the blades 
in the test section for the series. This is a prohibitive amount of work when the whole 
analysis is made using manual methods, but it may be possible to refine the technique 
and use machine ealculation to greatly simplify the work. 

The shedding of vortices from blades adjacent to the blade under consideration, 


is believed to have an effeet, of unknown magnitude, on the entire flow field, and 


14 


particularly on the blade considered. This would all tend to indicate the desired ana- 
lysis is of the whole flow field, the center blade in particular, with cross-plots of 
circulation to determine the effects of other vortex shedding on the circulation around 
this blade. 

In observing pictures of the sequence chosen, it will be seen that a considera- 
ble boundary layer exists along the last half inch of the nozzle blade. Also on some 
of the pictures a fairly large wake off the trailing edge of the nozzle may be observed. 
Both of these factors make it difficult to determine the exact edge of the blade itself. 
Although this edge was used as the zero point for measuring peripheral distance, a 
much better starting point would be one with a specific reference to a stationary object 
in the test section, notably, one of the vertical reference wires, or the crack in the 
glass of the test section, which runs horizontally between the blades. (As was men- 
tioned in Ref. 1, this crack does not affect the flow pattern, as can be seen by noting 
the continuation of the fringes, on both sides of it.) 

The centers of the fringes were established as accurately as possible, since 
these center points correspond most closely with the actual density, hence velocity on 
the blade. The boundary layer in the vicinity of the trailing edge presented some 
problems, but use was made of the assumption that the pressure is constant through 
the boundary layer, in a direction normal to the surface, to establish location of the 
fringes. This assumption is one conclusion of the order of magnitude analysis of the 
Navier-Stokes equation. 

The following assumptions were made in reducing the data obtained from the 


photo sequence, and are revicwed here: 


1.) Ff = fio = constant. Actually the highest velocity on the blade was 432 
ft./sec., which corresponds to the 6 1/2 fringe shifts. This is obtained from a value 
of F/e. = .933245, which indicates the degrec of validity of this assumption. 

2.) The flow is isentropic and the presence of the boundary layer is neglected. 

3.) All functions including V,?, and se arc continuous functions. 

4.) In Bernoulli's Equation, the term Pie for a givcn point is considercd a 
constant at any one instant in time. This makes it possible to find the velocity, using 
an iterative process. 

5.) The stagnation point on the leading edge of the blade is constant at a peri- 
pheral distance from the trailing edge of 2.73 inches. 

This last assumption was the object of an extensive analysis. Theoretically, 
the basis for finding thc potential function lies in the fact that one starts at a constant 
potential, preferably at an infinite distance upstream. However, due to the impracti- 
cality of this, the stagnation point on the leading cdgc was chosen. Since the location 
of this stagnation point was quite critical a lengthly analysis was made, thc results of 
which are tabulatcd in Table IX, and discussed in detail in the analysis scction of this 
report. Since the location of this point varied very little, it seemed logical to assume 
that the location was constant. This, in turn, besides simplifying the analysis, seemcd 
to reduce the fairly large variations in the graphical integrations which followed. On 
the other hand, in future analyses this assumption and the actual location of the stag- 
nation point should be chccked very thoroughly. Errors may actually be introduced 
by such an assumption. 


In addition to the above, an attempt was made to determine thc location of the 
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trailing edge stagnation point, which is a great aid in determining the actual velocity 
profile around the blade. From examination of the photographs little could be learned, 
except that the stagnation point is nearly constant, and located very near the trailing 
edge of the blade. The best indication is that the point is appEOMirentely PUZ inenes 
from the trailing edge on the pressure side of the blade. 

At this point in the analysis it was possible to check all the previous assump- 
tions, calculations, and results by comparing the velocity profile determined from 
analyzing the pictures with theoretical results. 

In Ref. 8, Bueckner and Schnackel, outline a method whereby the theoretical 
velocity profile around a turbine blade can be found with the aid of high speed compu- 
ters. The General Electric Company, used this method to find the theoretical velocity 
profile for the blades under consideration, with unity velocity upstream. The results 
of this program are plotted in Fig. 26 for an upstream velocity of 80.88 ft./sec., the 
same value determined for the picture series analyzed. Fig. 25 is a plot of V vs. S 
for one of the pictures, analyzed to show the type of results obtained in the analysis. 
Also plotted in Fig. 26 are several of the velocity profiles determined from other of 
the pictures analyzed. It is very interesting to note the excellent agreement in shape 
of the profile and in the actual velocities determined. If the theoretical curve is to be 
the guide, the initial assumptions and analysis are valid and yield excellent results. 

For reasons mentioned in the analysis section, only the suction side of the 
blade was considered in the integration involved to find oe - The method of graphi- 
cal integration used in the analysis was strictly a manual one. Since the velocity 


curve was continuous with no radical changes in slope, the mean ordinate for relatively 
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small values of s was multiplied by ds, beginning at S - 2.73 in. and progressing at 
even intervals toS = .25 in. It is also possible to use a planimeter or integraph for 
this type of work. Since the problem becomes one of iteration the best method, by 
far, would be to utilize a high speed computer. 

In the iteration process, values of Oe obtained as slopes from Figs, 44 - 53 
and Figs. 62 - 71,were taken considcring @ to he a positive value. Ref. 5 indicates 
this to be a valid conclusion. 

In the interest of simplification, and to reduce calculations, only cight values 
of V and oF /y¢ were selected from the curves in the first iteration. These corresponded 
to reference times, .16, .22, .28, .5-1, .40, .46, .52, .58 milliseconds. In the 
second iteration the intermediate times were used to give finite values of OPM many 
of which would have becn zero if the same times as for the first iteration had been 
used. Hence, the times, .19ms., .25ms., and .3lms. ctc. were used. It was be- 
lieved that this choice of different times for the second iteration may have induced 
errors in the calculations by shifting the velocity curve. To check this hypothcsis, 
the second iteration was recalculated using the samc times as for the first iteration, 
moo. Loms., .22ms., .28ms., re It was encouraging to note that these additional 
valucs did not change the curve, but more exactly defined the shape of the curve. 
This accounts for the fact that velocities were calculated for fifteen values of time 
in the second iteration. 

Figs. 72 to 74, are plots of V versus t for three values of 5, after two itcra- 
tions, The velocity for the zeroth iteration can be seen to vary for S = 1.90 inches, 


approximatcly 10 ft. /sec. which corresponds to a i/4 fringe-shift. For S = 2.50 
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inches, the velocity variation approaches 1/10 fringe shift. This accuracy is consid- 
ered to be quite adequate for an interferometric type analysis. However, as can also 
be seen in these figures, the variation for the first and second iteration appears to 
increase in amplitude. This divergence, and the factor of time, caused the analysts 

to discontinue the iteration at this point. However, the third or fourth iteration might 
have begun a converging trend. It is believed that more points for the same time 

span are necessary, to more accurately portray the true unsteady bchavior of the flow. 
If possible, it is also felt a more constant time interval between each photo is necessary 
However, this particular study was more of a feasibility experiment, and thus, the 
time variation between points is understandable. A possible way to meet both of these 
requirements is with a very high-speed Fastax-type motion picture camera which 
would give up to 12 or 16 frames per bucket passing sequence. The advantage of such 
a camera is that a true time study is made of the flow, where the assumption that the 
flow exactly repeats itself need not be made. 

At this point the photos were checked to note any correlation with the velocity 
plots in Figs. 72 to 74. The curves would indicate that a change in the flow field 
should be observed between t - .34 ms. andt =.37 ms. These times correspond to 
Figs. 10 and 11 respectively. Unfortunately no appreciable change could be noted, 
and no indication of any vortex shedding was evident. However, it is thought that the 
total blade must be considered, and more convergent velocities determined, before 
any comparison should be made with the interferomctric pictures. 

When the analysts could not explain the apparent divergence in the velocities, 


and the poor correlation with interferometric photos, a conference was held with Mr. 
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Hans Kraft and Mr. R. Fowler of Genera] Electric Company, in July 1960. At 
this conference,possible sources of errors were discussed. 

In the analysis, for reasons previously stated, the nose stagnation point was 
considered fixed in location on the blade periphery. It was noted this may result in 
small errors in the potential term, if a small area is neglected in the integration to 
determine the potential. However, it was assumed that this small error would be 
negligible when compared with the overall magnitude of the potential. But since ¢ is 
plotted versus time for a fixed peripheral distance, small plus or minus errors in the 
integration, due to shift of the stagnation point, could cause substantial changes in the 
@ versus t curve. This change in shape varies the slope of the curves and changes the 
magnitude of the Y?/st term. This slope is quite critical, since small changes in 
slope, for the small time intervals considered, may result in large variation in the 
magnitude of °*/¢. 

Due to the time element, the magnitude of this error could not be investigated 
further, but the data already obtained was give a cursory examination. The area in- 
volved is quite small, as the shift of the stagnation point was at the most, .05 inches 
on either side of the assumed stagnation point, 5 - 2.73 inches. The resulting error, 
however, may be compounded in the iterative process, and account for some of the 
large divergence in velocity from successive iterations. 

In order to determine the exact magnitude of any error induced it would be 
necessary to determine the exact location of the leading edge and trailing edge stag- 
nation points at all times. The approximate location may be obtained from the inter- 


ferograms if a fringe is located on the leading edge, but this would necessitate a 
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refinement in interferometric technique and a whole new series of pictures. This, 
plus the fact, that the techniques used in the analysis can not give the kind of accuracy 
needed to pinpoint the location exactly, indicate some of the problems involved. 

General Electric Company has recently determined the theoretical location of 
the stagnation point, for various flow conditions. Should a future analysis be made, 
the location of the stagnation point should certainly be considered, and experimental 
results compared with the theoretical data of the General Electric Company. 

If the iteration process would have converged to a final velocity, the analysis 
would have required the following steps. First it would be necessary to repeat the 
process for the entire blade, both suction and pressure sides. Having the final vel- 
ocities for the blade, the integral of V ds around the entire blade, will give the cir- 
culation around the blade. By noting the change in circulation with time, and the 
amount of that change, it may be possible to correlate the shedding of the vortices 


with changes in circulation, and to determine the strength of the vortices shed. 





CONCLUSIONS 


Several changes should be made in interferometric technique and in the analysis 

of the data before Ae results can be obtained. 

With regard to technique the following recommendations are made: 

(1+ More pictures per bucket passing sequence, should be obtained. 

(2) Equal time spacing, for each of the pictures in the bucket passing sequence, 
should be obtained. 

With regard to analysis of the data, the following recommendations are made: 

(1) The blade to be analyzed should be the center blade in the picture. 

(2) Methods of measuring the distance around the nozzle blade need to be 
simplified, and at the same time remain accurate. A template to the ex- 
act scale of the blade would save considerable time. 

(3) A complete analysis of the location of the leading edge and trailing edge 
stagnation points, and the variation with flow conditions and time, should 
be made. 

(4) Since the number of iterations required, and the type of analysis involved 
indicate a computer solution, the use of a computer to speed the analysis, 
and to more easily investigate the effect of variations in the parameters 
involved, is recommended. 

In general, three other recommendations are also made, 

(1) Fora more complete analysis, the flow fields around the other blades, 


and their effect on the blade being analyzed should be investigated. 
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(2) An investigation into the assumption that the flow repeats itself, should be 

conducted before future analyses are made with regard to vortex shedding. 

(3) An order of error analysis should be made to determine the magnitude of 

the error induced by neglecting the shift of the stagnation point. Such an- 
alysesshould also be made for errors induced in analyzing the interfero- 
metric picture, plotting the data and determining It) 54, 

Although as stated in the discussion, it was impossible to correlate the results 
of this report, with the interferometric pictures, this type of analysis appears to be 
feasible. Good correlation was obtained between the experimental and theoretical 
velocity distributions around the blade. The assumptions made appear to be valid, 
except for that of a constant leading Base stagnation point. Therefore, using the 
above recommendations, it should be possible to ultimately determine the actual flow 
around a nozzle blade, and to also quantitatively establish the strength of the shed 


vortex. 
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TABLE GH 


SAMPLE CONVERSION OF ROUGH DATA TO 
CORRECTED DATA 


on Photo Number 435 
Ref. Length .41 Inches 


Right Hand Blade Analyzed 














AS’ EAS” EAS Non Dimensional 
Between Points From Tip Corrected V V | 
Uncorrected Uncorrected & fSx aA Feet/Second V/80. 88 | 

06 04 359.3 3. 82 | 
. 41 ei 3960.58 Ae 
F tate) 64 414.7 4.41 
1.55 T.00s3 396.5 428 
1.76 L214 Sip ae. So2 
136 Lee oliso 3.38 
30 Zao 1.47 Age Oe!) Zag 
OU Zane 1.79 Zhows Loe 
0. 04 Fo 930) 140.25 1.56 
4.68 3.04 -140.29 -1. 56 
ayes, 3. 0U -159.0 -1.7305 
0. 66 3.68 -140.25 -1.,56 
6.86 4.45 -181.5 -1.98 
(ie es: 4.64 -215.3 -2.32 
71.44 4.82 -271.0 -2.99 
(Ey 4.91 =o1729 -3.38 
7.68 4.98 oo) ons -3. 82 
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TABLE III 


VELOCITY VERSUS PERIPHERAL DISTANCE FOR 


UM 435 

V 
. 04 359. 
5a 396. 
64 414, 
005 396 
14 Sag 
21 Shu 
47 2k 
cy ui 
30 140 
ie 0 
0-4 -140 
30 -159 
68 -140 
45 -181 
64 -215 
82 -271 
91 -317 
98 SSE 

t=-0ms 


Co Cc “1 Cf W 


ISO O& © © 


on 


PICTURE SEQUENCE UM 435 TO UM 458 


UM 436 
oe V 
. 06 Bio. oO 
ae 396.8 
OLS 414.7 
. 807 396.8 
911 359.3 
12390 ote g 
1.488 7a a oe 8, 
178350 215.0 
Does 140.25 
2213 0 
500 -140.25 
oS Ge -181.50 
sey | -140.25 
4.02 -159. 0 
2 es -181.5 
4.65 =2 0 Dies 
4. 84 -271.0 
4,93 -317.9 
4.98 -359.3 
t = .0707ms 
5S - inches 
V - ft/sec 


t 


- Relative Time, Milliseconds 
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UM 437 

S V 
ue 378. 
ae 336) 
eye, 414. 
.03 396. 
oe? OO. 
3 OL, 
ol pags: 
89 215 
42 140 
73 0 
94 -140 
34 oe 
fl -140 
Zs = ilo): 
60 S210 
89 Par): 
i) Sauer’ 
98 =309 
t = .0774ms 











Table II (Continued) 
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UM 438 
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378. 
396. 
414. 
396. 
359. 
oe 
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2D, 
140. 
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-140. 
-181. 
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-215. 
=z2il. 
-317. 
od), 


t = .1148ms 
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UM 439 

V 

.10 378 

“26 396 

67 414 

03 396 
16 Boos 
29 Bt. 
20 Zid 
87 2153 
34 140. ; 

ts Q 
93 -140. 
ao -159. 
62 -140. 
86 -159. 
40 -181. 
64 -215. 
83 -271. 
92 -317. 
98 -359. 

t = .1585ms 
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UM 440 

5 V 

10 378. 
ah 396. 
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wooo 396. 
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. 84 Zl. 
230 140. 
ee 0 
02 -140. 
28 -159. 
.61 -140. 
89 =) ele 
~A7 -181. 
arg | -215. 
. 86 -271. 
94 -317. 
. 98 -359. 
t = ..1753ms 
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t= .2159ms 
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70 ~215. 
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t = .2832ms 


UW OO WMOHwWH WON WM 


OoOomooenwn on ~ 
On 


UM 443 

5 V 

13 DiS. 

34 Sele. 

alee, 414, 
eo Sse 
Leo Saye! 
| se) S1F.. 
1.47 PA alee 
Loe VAISS 
220 140. 
2a te 0 
J 02 ~ 80. 
eo ee) - 80. 
4,00 - 80. 
4521 ~140. 
4.69 ~215. 
4.86 -271. 
4.94 -317. 
4.98 =O 

t= .38102ms 


wowowonwoanmd ms 


wo © OWN As =I] 
Oo 





213 378.5 .14 378. 
233 396. 


Table II (Continued) 
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UM 450 

S V 
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20 414.7 
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1.05 396. 8 
17 359.3 
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1.46 TaN) 
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2.30 140, 25 
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Table III (Continued) 
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5 V 5 V 5 V | 
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ie 2 317.9 16 359.3 £03 396.8 
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Ze i3 0 26 140.25 76 Zhou 
2.05 -140.25 
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4.47 -181.5 2 -~110. 
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Table II (Continued) 
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4.68 SA 1s Fas 
4.85 -~271.0 
4.93 -317.9 
4.98 -359.3 
t = .7669ms 
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TABLE IV 


CONVERSION OF REFERENCE DISTANCES 
TO REFERENCE TIME 


ee 











Reference Reference Reference Reference 
UM Distance Distance Distance Time 

Number (Uncorrected) (Corrected) From Zero Milliseconds 

Inches Inches Inches 






. 41 . 266 v 0 


20 130 PLoG 070662 
18 Phy 149 077417 
07 045 221 114827 
06 Ee 8 305 .158471 
vit ~ .0714 3375 .175306 | 
23 ~ ,1495 .4155 215885 
.43 2 (97g 5450 283171 
51 = 83l 597 310189 
60 = 4390 656 . 340844 
.68 - ,442 708 367862 
oa 550 . 166 397998 
.98 eG 902 . 468661 
01 - .655 921 478533 
214 - ,740 1. 006 . 522697 | 
19 ap 1.038 539323 | 
aa - ,850 1.116 , 579851 | 
35 =e 1.142 593360 | 
48 - ,96 1. 226 637004 
62 -1,05 ote 683767 
77 ~1,15 416 .735725 
. 87 SPI 1.476 . 766899 
01 -1, 305 1.571 . 816259 
. 06 =1,34 jh mts 834445 
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TABLE V 


VELOCITY AND POTENTIAL VERSUS TIME 
FOR SELECTED PERIPHERAL DISTANCES 


FIRST ITERATION 


UM § = 2.50 S = 2.20 > = 2.00 
Pictures Ref. Time V V V 
438 . 1148 126 169 LSS. 
439 wlo85 118 152 189 
440 ide 119 153 188 
441 noLoo 111 148 187 
442 POO! 108 138 173 
443 OLOZ 112 145 180 
444 . 3408 113 145 181 
445 . 3679 118 150 18s 
446 . 3980 113 144 1 
447 . 4687 109 136 175 
448 ~4785 110 141 182 
449 NAPA 118 ) sys) L933 
450 waooo PLS 148 184 
451 5799 120 LS: 191 
452 . 9934 109 141 178 
453 .6370 114 150 186 
5 -1.75 8 = 1.50 5 - 1.25 
438 . 1148 2o0 280 330 
439 E535 228 oes Bou 
440 ~1753 226 270 319 
441 WDA) 220 269 gaz 
442 p2000 282 256 320 
443 . 3102 218 Zoo ool 
444 . 3408 Zee 268 334 
445 Pood? 220 264 329 
446 . 3980 21S 262 335 
447 . 4687 22 21 336 
448 .4785 220 276 330 
449 Ds Ara 234 Zoe Bil F) 
450 2O390 224 269 Ses | 
451 5799 228 267 329 
452 . 59934 216 258 aio 


453 . 6370 223 264 321 








Table V (Continued) 








UM S = 1.00 eee 6 8 S =.50 
Pretures Ref. Time V V 
438 .1148 400 415 405 
439 was ogo 414 412 
440 Sires 395 405 405 
441 2109 398 407 409 
442 5 eae B90 415 408 
443 3102 398 415 408 
444 . 3408 404 415 410 
445 solves. Bo 414 411 
446 . 3980 403 415 408 
447 . 4687 403 414 411 
448 .4785 408 A 416 
449 Posed 403 425 424 
450 PoUO 403 420 420 
451 oOo 405 429 425 
452 Sac! oud 414 414 
453 Oo 402 434 427 
OE AO) 

438 . 1148 388 

439 ~1585 396 

440 EOS 395 

441 PoEOo 398 

442 P2ooc 388 

443 -o LOZ 390 

444 . 3408 B92 

445 . 3679 oot 

446 Sogou 390 

447 4687 390 

448 .4785 403 

449 Sa | 411 

450 .5393 408 

451 ono 410 

452 Be: 405 


453 .6370 411 





TABLE V-a 


VELOCITY AND POTENTIAL VERSUS TIME 
FOR SELECTED PERIPHERAL DISTANCES 


Biko mT EE RATION 


UM 5 = 2.50 5 = 2520 S = 2.00 
Pictures Ref. Time D DP QD 

438 1148 17.54 Donte 100. 04 
439 joes: 16.45 50.45 92.95 
440 Ploo 16.65 50. 90 93.40 
441 e2o9 ba 92 48,42 90.17 
442 e2o02 15. 35 45.60 84.10 
443 woOL U2 EG. 30 48.955 89. 30 
444 . 3408 16.61 48.11 88. 86 
445 nooo 17.14 50.64 92.14 
446 - 3980 16. 67 48.92 88. 92 
447 . 4687 1 Gs22 46.62 86.62 
448 24785 Je asi 47.58 88.85 
449 mOeea 16.62 HUS6Z 04 4 
450 oso 16n72 49.46 90. 96 
451 .5799 eee Nal Dies G 94,36 
452 Povo 16. 41 47.41 ol o)e pedtl 
453 .6370 1O7 02 DO be A ONT 

5 = Lato see AG, So leo 
438 1148 ESS aay: 217 aoe Maen Ue 
439 roo 144.95 2075 10 284.70 
440 aloo 145.15 206.65 219.900 
441 Dood 141.92 203,42 219.42 
442 neooe P32..156 290,35 261.60 
443 Pow? 139. 05 199.05 213.00 
444 . 3408 139. 36 200.61 275.61 
445 SEXES) 142. 64 202.89 278.64 
446 . 3980 ERs he Poa 270. 66 
447 . 4687 Loo. 12 196.47 At ia Wee Fi 
448 .4785 138. 58 201.08 Zhao 
449 oad LAy LZ Poe 288. 37 
450 4 abs 8 141.71 203.46 210. Ui 
451 ae, 146.61 208. 61 281. 86 
452 . 0934 136.16 195,16 260,06 


493 6370 143.37 204.37 21 G.62 





Table V-a (Continued) 


UM 
Pictures 


438 
439 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
450 
451 
452 
453 


438 
439 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
450 
451 
452 
453 


Ref. Time 


. 1148 
~ 1585 
too 
aloe 
-2002 
3102 
. 3408 
. 3679 
. 3980 
. 4687 
4785 
0220 
. 5393 
5799 
0934 
. 6370 


1148 
. 1585 
1753 
, 2159 
2 Oa7 
. 3102 
. 3408 
pallies 
. 3980 
4687 
4785 
.o2at 
.0393 
5099 
. 9934 
. 6370 





200 





TABLE VI (1) 


DETERMINATION OF CORRECTED 
VELOCITY BY FIRST ITERATION 


me rem arnt ee tn a i ee Se a te ee ee 


























t =-.16ms 
V v2 JO/d2 
2. 0 0 
2. 118.0 13,924 - 1,540 15, 464 124.3 
2. 152.7 23,340 = 5,660 29, 000 170.3 
2. 189.3 35, 830 ~ 9, 820 45,650 213.7 
‘. 22.9.0 52, 441 -13, 020 65,461 255.8 
ie 273.0 74,529 ~14, 450 98,979 314.7 
r, 323.2 104, 329 -16, 120 120, 449 347.0 
iu 397.8 158, 200 -24, 850 133, 050 427.9 
75 409.6 167,750 -22, 229 139, 972 435.7 
50 408. 6 166, 850 -22, 222 189, 072 434.8 
25 395.0 156, 025 -18, 050 174, 075 417.2 
22 
V ve Oy Vz OF Vcorr | 
2. 0 0 0 0 0 
D. i 3 12,400 - 625 13, 025 114.1 | 
2, 145.0 21, 025 = Ba soe 24,585 156.7 | 
2. 181.5 32,900 ="2s 360 40, 260 200.5 | 
1. 219.6 47,250 “li 0 58, 350 241.6 
i 263.5 69, 350 11400 80, 750 284.1 
i), 320.5 102,700 -14, 600 117,300 342.5 
1. 396.0 156, 816 -14, 600 Minh, 254g 414.0 
407.6 166,100 -14, 850 180, 950 425.4 
408, 2 166, 600 -16, 800 183, 400 428.3 
397.8 158, 200 -23, 300 181,500 426.0 





| 
| 
j 
| 


a ri er ae nn = 








TABLE VI) 















t- .28 
V v2 ra re 

Zs 0 0 0 
7 108.2 11,720 973 
2. 139.0 19,321 -1,500 
Ze oe 0 50,916 0 
i. es 2 45, 000 Soo 
i, 296.1 65, 600 -1, 805 
1. 3o20.0 106, 000 0 
i BIO. 156, 600 =12529 

~19 414.2 171, 400 0 

50 408. 0 166, 500 -7, 000 

Zo Soleo 150, 900 4,360 

V v2 

Za aS Q Q 
2.50 113.5 12, 880 1, 250 
2. 20 146.5 21,450 Oo, 450 
2.00 To2. 7 33, 400 7,980 
1.75 aS 49,700 10, 410 
£.50 201.0 iL. GOG 8,680 
LAS. Soles 110, 000 9,450 
1.00 403.5 162, 800 6,530 

19 415.0 Li 2.220 4,720 

20 410.0 168,100 6, 600 

20 a71l.8 153,500 7,640 


ee ee ee ee 


106, 000 
158, 129 
171, 400 
173, 500 


145,140 


62, 920 
100, 550 
156, 270 
167,505 
161, 500 


145, 860 





SOnnunn kt Oh Cl 


0 
6 
neo, 
ae 
a) 
<0 
3 
ar 
20 
io 











TABLE VI (3) 














Uae vee 
V Me O91 Oe 522% 

2. 0 0 0 0 
Ze 114.5 13,100 carl, 14,420 120.0 
2. 145.0 21,025 - 5,000 26,025 161.1 
2. 176.5 31,150 -11, 650 42, 800 206.9 
a. 214.8 46,100 e200 67, 800 260. 4 
Ly 261.8 68, 500 -14, 450 82, 950 288.0 
1, 334.6 111,900 -15, 300 127,200 356. 8 
Me 401.6 161, 250 - 7,900 168, 750 410.7 

414.1 171, 450 aM!) 177,410 421.2 

409.5 167,600 - 5,280 172, 880 415.7 

389.6 151, 800 -20, 900 172, 700 415.6 | 





2. 0 0 0 0 
2. 109.6 12,000 0 12, 000 109.6 | 
2. 140.8 19, 810 833 18, 977 MEY TE 
2 174.5 30, 420 3, 750 26, 670 163.5 
ie 218.0 47,524 3,190 44,344 210.6 
ip 265. 8 70, 650 5, 980 64, 670 254.2 
iif 335.6 112,600 16,120 96, 480 S107) 
I. 404.1 163, 300 5, 980 157, 320 396.8 | 
415.4 172, 500 8,330 164,170 2 
411.4 169, 300 10, 000 159, 300 ei 
389.5 4 


151,700 15,280 136, 420 


Oe ge ee | ee ee ey — pee - - ee ome - —.- —— ewe 





TABLE VI (4) 












t= .52 
V ve Ob/e 

as 0 0 0 
Ze 114.6 13,150 906 12, 594 ie ZL 
a 149.6 22, 400 4,100 18,300 130.2 
Ze 188.6 35, 600 2,910 32,690 $3.5 
A Zooaee 54, 600 3000 00,750 Z2IeoS 
i, 282.0 79, 524 1, 667 1, oot 279.0 
1. Sool 112, 350 7,500 104, 850 Aah (3) 
3 403.6 162, 900 4,580 158, 320 Byelt| ccs) 

ee a235.0 179,300 3,470 175, 830 419.2 

50 420.8 ~ 0 - he 

25 «409.5 a 0 = 409.5 

t- .58 
ID/4 UV oGAe Vceorr 

Belo 0 Q Q 0 0 
2.50 TET .3 - 0 - LAT 23 
eae 153.4 Zonas U 417 Zo, lis sy 7a 8, 
2.00 186.3 34, 700 ~5,140 39, 840 E9925 
He 15D A | 49,825 -5, 000 a4, 825 234.2 
1.50 261.2 68, 200 ~8, 890 (A CSUSAG, VAL AT ey 
25 326.4 106,500 ~9, 880 116, 380 341.2 
1.00 404.5 163,600 ~11,670 175, 270 een | | 

is 423.9 179, 500 -9,310 188, 810 434.7 

50 424.2 179, 900 -7,230 18454350 432.6 

29 408.4 166,600 -7, 920 174, 520 417.7 


are — oe er ee 





wa 
Lo Mle sec 


Ref. Time 


.16 
a2 
. 28 
34 
. 40 
» 46 
2 02 
. 98 


Ret. Time 


.16 
Sag 
P20 
. 34 
. 40 
. 46 
202 
209 


Ref. Time 


SLO 
22 
.28 
. 34 
» 42 
. 46 
. 92 
. 98 


VELOCITY AND POTENTIAL VERSUS TIME 
FOR SELECTED PERIPHERAL DISTANCES 


TABLE Vil 


SECOND ITERATION 


S = 2.50 S$ = 2.25 
AVS ons. \ ee 
124.3 17.81 170.3 54.81 
ieee a73 156.7 50.53 
103.5 15.16 144.2 46.91 
108.0 16.54 126.6 45.04 
1Gron 744 161.1 52.59 
109.6 16.42 137.7 46.92 
Moet 7.04 135.2 47.04 
hie7 2 tro 152.0 50.24 
§=1.75 S$ =1.50 
V D V ¢g 
255.8 161.81 314.7 232.56 
241.6 150.03 284.1 215.53 
215.4 135.41 260.6 194.41 
198.1 125.29 250.9 182.04 
260.4 155.09 288.0 223.84 
210.6 133.67 254.2 191.92 
229.5 138.19 279.0 202.29 
234.2 147.49 277.7 210.99 
S = 1.00 S=.75 
V gp V dD 
427.9 417.31 435.7 525.81 
414.0 389.78 425.4 495.03 
397.6 362.41 414.2 464.41 
395.3 342.92 409.3 444. 42 
410.7 401.84 421.2 506.34 
396.8 352.67 405.2 453.42 
397.8 367.66 419.2 470.66 
418.7 385.24 434.7 492.24 


S = 2.00 
ae 
213.7 £102.81 
200.5 95.03 
176.0 86.91 
TO9eo 80.29 
206.9 98. 09 
163.5 85. 67 
183.5 87.04 
£99.35 93.74 
5 =1.25 
22 
347.0 316.06 
342.5 293.53 
325.5 270.41 
317.0 252.54 
356.8 305.59 
310.7 262.42 
G2o.cue 21 to 
341.2 287.74 
5 = .90 
V DP 
434.8 634.81 
428.3 601.78 
416.6 568.66 
401.8 546.42 
415.7 610.59 
399.1 554.67 
420.8 576.16 
432.6 599.74 





TABLE VII (Continued) 


S = .25 
Ref. Time V fp 
.16 Ae ee, (41.56 
aoe 426.0 708.53 
. 28 381.0 670.41 
34 381.8 644, 67 
40 415.6 714. 34 
46 369.4 651.42 
52 409.5 680.16 
58 417.7 706.24 








TABLE VIII (1) 


DETERMINATION OF CORRECTED 
VeBoOcily BY SECOND ITERATION 


























ens | 
Va Vaz Oo). VV O- Oh ye Veorr | | 
0 0 " : = | 
118.0 13, 924 0 14,674 121.1 
152.7 23,340 - 7,500 30, 840 175.6 | 
139.3 35, 830 oe Teas 46, 930 216.7 
229.0 52,441 -13, 600 66, 041 257.0 | 
273.0 74, 529 -18, 200 92,729 304.5 | | 
323.2 104,329 ~21, 400 135,729 368. 501 | 
397.8 158, 200 -29, 900 138,100 433.7 | 
409.6 167, 750 -39, 200 206, 950 454.8 | 
408.6 166, 850 -51, 400 218, 250 467.2 
395.0 156, 025 -35, 100 191, 425 137.5 | 
Va Va" 
0 0 
all 12, 400 - 3,680 16, 080 126.7 
145.0 DUN 25 = 5a550 VAS aos 163.0 | 
181.5 32, 900 -11,100 44, 000 209.8 — 
219.6 47, 250 ~17, 350 64, 600 254 
263.5 69, 350 -29, 700 99, 050 314.7 4 
320.5 102,700 “317.00 134, 400 366.7 | 
396.0 156, 816 -41, 700 198, 516 445.5 | 
407.6 166, 100 -45, 000 211,100 459.5 | 
408. 2 166, 600 -46, 100 212,700 461.2 | 
397.8 158, 200 -52, 500 210, 700 459.0 | 











TABLE VII (2) 











re Oe 
V i IP/54 We OU e Vcorr 
2.73 0 0 ks : = 
2.50 108.2 11,720 - 570 12,290 110.9 | 
2.25 139.0 19, 321 = 3 BNW 22,931 151.3 | 
2.00 176.0 30, 976 ~11,100 42,076 205 ie 
L. 75 212.2 45, 000 -23, 100 68,109 260.9 | 
1.50 256.1 65,600 -24, 400 90, 000 300.0 
1,25 325.5 106, 000 -28, 500 134, 500 366. 8 
1.00 395.5 156, 600 -35, 200 191, 800 437.9 | 
75 414.2 Tyee oe ERA’ 210,100 458.5 | 
50 408.0 166, 500 -37, 500 20-4, 000 451.7 | 
25 387.9 150, 000 -~45, 000 195, 000 441.7 | 





2. 
2. 9 
2. 146. 21, 450 S505 21,955 148.2 
2. 182. 33, 400 - 6,250 39, 650 194.0 — 
1, 222. 49,700 3, 190 52, 890 230.0 | 
i 267. 71,600 0 71, 600 267.6 — 
i 331, 110, 000 0 110, 000 331.8 | 
L. 403. 162, 800 0 162, 800 403.5 | 
415. 172, 225 0 172, 225 415.0 — 
410. 168, 100 0 168, 100 410.0 
391 153, 500 0 is. ty 








TABLE VIII (3) 











t = .40 

v2 OF/ + 

0 a 

13, 100 0 
21, 025 4,170 
31, 150 2, 430 

46,100 Q 

68, 50 0 

Pit 800 0 

LGi, 250 0 

171, 450 0 

167, 600 " 

151, 300 0 

t= .46 

V v2 MYL 

2 0 0 - 
2.50 109.6 12, 000 - 500 
a. 20 140.8 19, Si0 - 2,710 
2.00 174.5 30, 420 ~ 8,470 
t.7D 218.0 47, 524 a2 O20 

1.50 269.8 70,650 0 

1.25 335.6 112, 600 0 

1.00 404.1 163, 300 0 

79 415.4 72.500 0 

20 411.4 169, 300 0 

5) 151, 700 0 








a a re NiCOLE 
13, L00 114.5 
16, 855 129.9 
28,720 169.5 
46,100 214.8 
68,500 261.8 

111, 900 334.6 

161, 250 401.6 | 

171, 450 ee 

- 409.5 
2 a 
ve Ib, | 
12, 500 111.8 
22,520 150.0 
38, 890 197.2 | 
50, 444 224.6 | 
70,650 265.8 

112, 600 335.6 

163, 300 404.1 | 

172, 500 415.4 | 

sz 411.4 | 
Bh, ele 








TABLE VU (4) 
























= 52 
ee ee oe eee 
2. 0 " : ; 
2. 114.6 13, 150 625 12,525 0 
2. 149.6 22, 400 2,640 19, 760 7 
2. 188.6 35, 600 7,640 27, 960 167.4 | 
i, 233.8 54,600 12, 200 42, 400 205.9 
1. 282.0 79,524 17, 100 62,424 249.9 : 
ie 335.2 112, 350 22, 500 89, 50 299.8 
i 403.6 162, 900 23, 300 139, 600 373.7 | 
15 423.5 179, 300 28,100 151,200 388.8 | 
50 420. 8 177, 000 32, 500 144, 500 380.1 | 
25 409.5 167,700 45, 000 122,700 350.3 | 
S056 
V y2 D y At < J. V 
VV Be corr 
2. 0 0 : - 
2. 117.3 13,760 139 13,521 116.2. 
2. 153.4 23, 530 6, 390 17,140 131.0 | 
2: 186.3 34, 700 11,100 23,600 yoy 
iV 223.3 49, 825 10, 830 38, 995 197.5 
i. 261.2 68, 200 4, 580 63, 620 252.3 
1. 326.4 106, 500 4, 450 102, 050 319.5 | 
1. 404.5 163, 600 25, 300 138, 300 alee 
423.9 179, 500 31, 700 147, 800 384.5 | 
424.2 179, 900 32, 500 147, 400 384.0. 
408.4 166, 600 31, 900 134, 700 367.0), 


nee ee — ee ae = = - son ee A a ee ee — ee 


| 
| 
| 
| 








TABLE VII (5) 

















V 
e203 0 0 0 0 0 | 
a 0 12, 996 - 1,250 14, 246 Noses: : 
peo 148.0 21, 904 - 5,980 27, 884 167.0 | 
2.00 185.4 34, 340 ~10, 820 45,160 22 <5 : 
me. 75 224.5 50, 350 ~16, 400 66, 750 258.4 | 
1,50 268.0 71, 824 -23, 900 OGna24 309.4 | 
125 320.0 102,400 ~31, 900 134,300 366.6 | 
1.00 396.8 157,400 ~38,100 195,500 442.2 | 
164, 500 ~41, 700 206, 200 454.1 
165, 200 -45, 800 211,000 459.4 
158,000 -45,100 203,100 450.6 
t=.25 ans ae 
ve OF V8  Veorr | 
E73 0 0 0 0 0 
2.50 109.4 11,950 -~ 2,710 14, 660 121.0 
2.25 et 5 19, 950 - 4,660 24,610 156.9 
2.00 177.8 31, 500 -10, 820 42, 320 205.8 
75 214.9 46,150 ~20, 600 66, 750 25302) 
1.50 259.5 67,300 -29, 200 96,500 310.7 | 
125 319.6 102,100 -30, 600 132,700 364, 4 | 
1,00 395.4 156, 200 -38, 900 195,100 441.7 | 
75 411.6 169, 500 -41, 700 211,200 459.6 
50 408, 2 166, 600 -43, 800 210, 400 458.8 
25 POUL S 153, 500 -51,400 | 204,900 452.6 








TABLE VIII (6) 








t -.3l1 ms “| 
a none | } a 
Sa : vm Oty), oe 27. " Veorr 
0 0 0 0 0 
110.0 12,100 2,500 9,600 98.0 
142.0 20,164 ~ 2,500 22,664 150.3 
177.0 31, 329 - 8,750 40, 079 200.1 | 
216.0 46, 656 =17,100 63, 756 252.4 | 
254.9 65, 000 -18, 600 83, 600 289.2 
B202D 107 200 -26, 400 133,600 BO0uS 
398.2 158, 500 -28, 350 186, 850 432.3 
415.0 172,225 =3 300 205, 525 453.2 
408.5 166, 800 ~34, 200 201, 000 448.3 | 
390.1 152,150 _ _-39,700 ———s:191,850 438.0 | 
Losrstot MS 
Vv v2 a iy 
0 0 0 
ley 5 13, 800 1,210 12, 590 Nitalgae | 
149.6 22, 400 10, 410 11, 990 109.5 | 
183.2 33, 550 24,700 8, 850 94.0 | 
222.0 49,284 41,500 7,784 88.3 | 
264.8 70, 000 58, 100 11, 900 109.4 | 
333.6 111, 300 73, 600 37,700 194.1 | 
402.8 162, 200 82, 000 80,200 283.2 | 
414.5 171, 800 86, 000 85,800 292.9 
410.4 168, 400 89, 200 79, 200 281.5 _ 
391.0 152,881 ss 96, B00 56,081 —«- 236.8 | 








TABLE VI (7) 


+) 


en no NOM NOME OE OE 


ee 


© 
© 


OW FR KY wn OK WO CO 


Oo Oo — © OO NM © oF fe © 


166, 490 


45, 900 





120, 550 


Vcorr ; 
0 0 0 0 
12, 321 - 2,220 14,540 120.5 | 
19, 750 - 7,780 pad (neo pil!) 166. 0 | 
30, 450 -17,100 47,550 218.0 | 
45, 700 -30, 000 fone Dae Z | 
68, 000 -44, 000 112, 000 334, 6 
112, 400 -59, 000 171, 400 414.0 | 
161,650 -67, 800 229, 450 479.0 
171, 450 -72, 300 243, 750 493.0 
167, 400 -73, 600 241,000 491.0 | 
1Si.300 -88, 900 240, 200 490.0 | 
49 ms =| 
v2 OD/ >, vee, Vcorr | 
0 0 0 0 : 
12, 321 820 11,501 107.5 
19, 700 0 19, 700 140.4 
33, 600 2.220 31, 380 bid. ul 
52, 800 13, 900 38, 900 197.3 
78, 900 14, 600 63, 900 252.8 
112,800 19, 900 92, 900 304.8 
164, 000 22,200 141, 800 376.6 
i430, £00 23,600 154, 500 393.0, 
173.000 32,600 151,300 389.0 








TABLE VII (8) 


167,300 


128,400 


7 | t= 55 ms i. a 

S V v2 Od; » 4 V cm OU 6 Veorr | 
2.73 0 0 0 0 0 : 
9.50 117.0 13,700 347 13, 353 ils Ge 
2.25 154.6 23, 850 4,300 19,550 139.9 | 
2,00 188.6 35, 500 9, 450 26, 050 161.2 | 
1,75 231.0 53, 360 12, 630 40,730 201/38 
1,50 271.8 73, 800 12,500 61, 300 237 6a) 
1.25 332.8 110, 700 15, 280 95,420 309.0 | 
1.00 403.3 162,600 23,900 138,700 372.5 | 
75 421.6 177, 600 29, 600 148, 000 384.7 | 
50 423.0 178, 900 30, 900 148, 000 S340 | 
25 409, 0 | 38, 900 358. 4 | 
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ig Sal Sn GOL‘S FEV Sak 89°S SSP 

SO ae GL°S CEP SHA GL°S LSP 

00 °T- CFL°S o&P SHA 9L°S 9SP 
096° - CSL°G TEP SHA CLS CSP 
8E8° - GL°SG OFF PLS PSP 
oo = De eG 660 G89°SC ESP 
ton = GL°S SGP SAA ee < CSP 
¢69° - VLG LGF SHA 9L°S TSP 
Coa = PLS 9GP SdA LL°S OSP 
OLS” - PL °S GCP SAHA CL°S 6PP 
O9T* - ta 6 PCP SHA GCL °C 8hP 
Ort = CL°S ESP SHA PLS LEP 
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86G° - 91°C TOP SAA STL°S CVT 
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Figure 12 
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Fizure 17 
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Flg. 57 velocity versus peripheral pistance for t = .54 milliseg, 
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Fig. 58 velocity versus peripheral pistance for t = .40 millisec, 
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Fig. 61 velocity versus peripheral pistance for t = .58 milligec 
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Fig. 65 velocity and potential versus Time for gs — 1.00 inches 
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Fig. 71 velocity and potential versus pime for g = 2.50 inches 
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APPENDIX 


For the benefit of those readers who do not have access to Ref. 1 and 2 the 
following excerpts are included to give a better understanding of the principles of in- 
terferometery, the test set-up involved, and the method of data reduction. 

INTERFEROMETERY 

Fundamentally, the optical method used in the interferometer is based on the 
fact the speed of light depends on the index of refraction of the medium through which 
it passes, and the index of refraction of a gas in turn depends on its density. Ina 
high-speed gas flow the density changes are sufficiently large to make these phenomena 
sizable enough for optical observation. The interferometer measures directly changes 
in density, and is primarily suited for quantitative determination of the density field. 

Fig. 1 shows the essential parts of the Mach-Zehnder interferometer, the type 
used in this investigation. Referring to Fig. 1, light from the source first passes 
through a collimating lens which renders the light parallel, and then passes through a 
monochromatic filter. It then passes through the first "splitting plate" which is a 
plane half-silvered mirror that passes half the light and reflects half the light. The 
light which is reflected is changed back to its original direction by the upper mirror 
and then passes through the test section. The light which passes through the first 
splitting plate passes through a pair of glass compensating plates which duplicate in 
thickness and quality the glass side walls of the test section, after being directed to 
the upper splitting plate by means of the lower mirror. The two beams which divided 


at the lower splitting plate are recombined at the upper splitting plate and are then 





focused by a lens system on a photographic plate. The chicf difference between the two 
beams of light is that the upper beam has passed through the test section. Since the 
different rays of the upper beam are retarded by different amounts as they pass through 
portions of the test section of different density, the rays will have various phase dif- 
ferences with the corresponding rays of the lower beam when they are recombined, 

and thus an interference pattern will be formed at the photographic plate. 

The interference pattern is composed of contours of brightness. These con- 
tours represent contours of constant index of refraction, and therefore represent 
contours of constant density in the flow field. Fig. 1, of the report, shows such an 

"infinite-fringe" interferogram of the flow through the turbine nozzle, the term "in- 

| finite-fringe'' signifying that the light field is uniform in the absence of flow through 
the test section. The dark bands in the figure are loci of points where there is com- 
plete interference, and represent contours of constant density. 

| For a more accurate quantitative evaluation, the fringe-displacement method 

1 may be used. If the second splitting plate is rotated through a small angle with respect 
to the first splitting plate, two ''coherent" beams of light which were in phase at the 
first splitter will, through the change in lengths of light paths, be out of phase at the 
screen. If the splitters are rotated about axes normal to the beam direction, there 
will be formed on the photographic plate, or screen, successive light and dark fringes, 
uniformly spaced, with cach fringe lying parallel to the axis of rotation. 

If the air density in the test section is uniformly increased, this will produce 
a uniform displacement of all the wave fronts passing through the test section. This 


displacement in turn will cause the interference bands on the screen to shift in a 











direction normal to the bands, although the bands will remain parallel and uniformly 
spaced. The fringe shift is a measure of the change in density in the test section. 
When there is flow in the test section, with corresponding nonuniform density changes, 
similar fringe shifts will occur, but, as they are no longer uniform, the resultant 
fringes will be curved. 

To analyze the interferometric pictures it is necessary to determine the den- 
Sity associated with each of the fringes, since changes in density are all the interfero- 
meter can measure. 


The free stream density § is determined from the values of Pg, T. 


and cis 
Which are either determined during the test runs or calculated from the test data. 


The density of the remainder of the flow field is determined from the two-dimensional 


fringe displacement relation: 


= 1+ X (Ae) Fo 





Where; 
9 - local density of the air (lb/ft) 
Se = Reference density of the air in the undisturbed fringe field (1b/ft?) 
y = Fringe displacement normal to the undisturbed fringe direction (in. ) 
b = Fringe spacing (in) 
A. - Wave length of light used, in a vacuum (mm) 
os = Length of light path in turbine test section (mm) 
net Gladstone-Dale-Constant for air 


a 8 
i 


This is the relation used for the reduction of data from the fringe displacement 


pictures. The data reduction of this type of interferogram is very tedious, time 





consuming, and prone to error, while the use of the infinite fringe interferogram 
yields constant density contours on the basic interferogram. The only data reduction 
that is required is the determination of the change in density between any two isopycnic 
(constant density) contours (which is a constant for the whole field if the properties of 
the interferometer optics are precise enough) and a calculation of the density in the 
undisturbed flow region just upstream of the nozzle blade row. Fortunately this den- 
sity, Se has been found to be nearly constant with time, which is desirable for com- 
parison purposes. 

For the series of pictures considered in this report, the density can be deter- 
mined by the following relation obtained from interferometric theory and utilizing 


light of a given wave length (4358A°) with a specified path length of light in the turbine 


af - 001027 


The above relation gives the ratio of the change in density, A Sreferred to the density 


test section: 


of the undisturbed upstream field, Se between any two isopycnals. 

If each dark fringe is numbered in sequence starting with number one at the 
inlet to the nozzle blade row, the density ratio SY 9 ot each isopycnal can then be 
determined. 

It should be noted that since the upstream field is light the density ratio of the 


_ 9. 01027 while the AG 9 between any two 


first isopycnal would be ( S/ ah = | 
succeeding isopycnals would be equal to 0.01027. This reduction of data is well 


illustrated in Tables I and II. 








TEST SET-UP 

The radial flow, single stage impulse turbine used for this study was designed 
and fabricated by the General Electric Company. A cross sectional and side view of 
the turbine appear in Figs. 3a and 3b respectively. 

The turbine Wheel has 208 bucket blades mounted in the periphery. The outside 
diameter of the wheel at the buckct leading edge is 41.75 in. The span of the bucket 
blades is 6.00 inches. Four cutouts spaced 90° apart are located in the end plates 
that contain the buckets to make visible the air flow in these four bucket channels. 

The bucket cutouts are terminated by one-inch diameter holes downstream of the 
moving stage to allow an opening for the beam of light used by the clectronic flash 
triggering apparatus. The bucket section is General Electric Type T~-1984237 with a 
27° bucket entrance angle and a 24° bucket exit angle. With the buckets evenly spaced 
along the periphcry of the turbine wheel the bucket pitch spacing at the O.D. of the 
Wheel is 0.631 in. The pitch spacing at the trailing edge is 0.600 in. 

The partial admission air inlet section consists of a wooden bell mouth faired 
into a stcel rectangular flow passage (6.35 in. x 18.98 in.). Provision is made for 
mounting twelve nozzle partitions with six of the nozzle partitions mounted between 
optically flat glass plates for viewing purposes. These six steel nozzle partitions 
are isolated from the glass plates by plastic gaskct material to eliminate local stresses 
in the glass. The nozzle partition is General Electric Type K-6915600. The nozzle 
sections are located on a basis of 104 equally spaced in 360° with the nozzle trailing 
edges locatcd on a 41.880 in. diameter giving a nozzle pitch spacing at the trailing 


edges of 1.295 in. The nozzlc exit angle is 13°59) he partial admission assembly can 





be raised or lowered with respeet to the turbine casing giving a variable gap between 
the nozzle trailing edges and the bucket leading edges. 

The air that passes through the turbine is inducted from the laboratory as a 
stagnation reservoir. The air is taken from a height of about ten feet above the lab- 
oratory floor through a calming section which is connected to the turbine partial ad- 
mission inlet seetion. 

The air flow then passes through the turbine stage and is turned and passes 
through the turbine outlet in an axial direction. 

A sketeh of the air induction system with the turbine is shown in Fig. 5. In 
the Figure a blower is shown at the turbine outlet and was used for inducting air into 
the turbine. The blower has been replaced by an air ejector arrangement which has 
increased the flow rate through the turbine, with two attendant advantages. First the 
increased density change through the cascade caused a much greater fringe shift on 
the interferogram, yielding more accurate quantitative information, and second, 
placed a sonic throat between the driver and turbine section eliminating any unsteady 
effects propagated upstream by the air induction system. 

The turbine load consists of a10 H. P., D.C. generator and a bank of 25, 300 
watt lamps. The load is adjusted to the desired value by means of a variable rheostat 
in the generator field and by selecting any desired number of lamps. 

ANALYSIS 

For general information Figs. 2 and 4 are included. Fig. 2 is a sketch of the 

vortex shedding phenomenon, showing the general location of the shedding and the 


path after shedding. UM Figs. 435, 451, 453 show some good examples of this vortex. 











Fig. 4 illustrates the numbering of the fringes and the direction of measure- 
ment of the peripheral distance. From the numbering can be gained an idea of which 
fringes correspond, and the location of the fringe of highest density. The right hand 


blade shows the location of certain peripheral distances mentioned frequently in this 


report. 
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